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Report  Summary 


In  this  project,  we  have  developed  a  compliant  mechanism  synthesis  theory  that  incorporate  a  gen¬ 
eral  framework  for  determining  pseudo-rigid-body  models,  type  synthesis  algorithms  (determining 
mechanism  topology  for  a  specific  task),  kinetostatic  analysis/synthesis.  This  research  starts  from 
classification  and  type  synthesis  of  flapping  wing  MAVs.  We  have  studied  over  15  flapping  mech¬ 
anisms  in  the  literature  and  classify  them  based  on  workspace  type,  linkage  topology,  rigid  body 
or  compliant,  actuator  type  and  mobility.  This  classification  lays  out  the  foundation  for  type  syn¬ 
thesis  of  flapping  wing  mechanisms.  We  have  also  developed  a  parameter  optimization  framework 
for  determining  the  pseudo-rigid-body  (PRB)  model  of  cantilever  beams.  A  novel  concept  of  PRB 
matrix  has  been  proposed  to  describe  the  general  topology  of  an  arbitrary  PRB  model. A  novel 
3-spring  pseudo-rigid-body  model  for  soft  joints  with  significant  extension  effects  was  proposed. 
These  soft  joints  are  often  used  in  transmission  mechanisms  in  flapping  wing  micro- air  vehicles. 
Based  upon  this  synthesis  theory,  we  have  shown  that  compliant  joints  can  significantly  reduce 
the  power  consumption  in  flapping  wing  MAVs  if  they  are  designed  to  enhance  the  down  stroke 
flapping  motion.  We  have  implemented  this  theory  into  a  computer-aided  design  software  called 
DAS-2D  with  a  rich  user-interface.  It  has  been  applied  to  solve  several  structure  design  problems 
including  design  optimization  of  compliant  transmission  mechanism  for  flapping  wing  MAVs,  a 
bistable  buckling  beam  design  and  a  robotic  finger  actuated  using  shape  memory  alloy  actuators. 

The  accomplishments  and  new  findings  of  this  project  are  listed  in  below: 

1.  We  have  conducted  a  comprehensive  survey  of  existing  MAVs  from  around  the  world.  In 
particular,  we  have  categorized  about  15  MAVs  based  on  workspace,  compliant  or  rigid  body, 
type  synthesis,  mobility,  and  actuator  type.  This  survey  is  expected  to  serve  as  a  resource 
for  continued  design  and  development  of  smaller  and  more  efficient  MAVs. 

2.  We  have  built  a  comprehensive  dynamic  mathematical  model  of  the  entire  drive  mechanism. 
To  this  end,  Theodoresens  theory  for  predicting  the  lifting  force  on  an  oscillating  airfoil  was 
chosen  to  calculate  the  lift  force  on  the  wing.  A  sensitivity  analysis  was  performed  on  the 
model  to  study  the  effect  of  system  variables  on  the  performance  parameters:  lift  produced, 
power  consumed  and  lift/power  ratio.  The  critical  parameters  were  identified  and  run  through 
an  optimization  algorithm  to  explore  the  design  space  for  improving  the  performance  of  the 
MAV.  It  was  seen  that  the  power  consumption  was  reduced  by  73.8%  without  affecting  the  lift 
significantly.  The  results  also  suggest  that  some  parameters  affect  lift  and  power  differently, 
and  in  general,  the  variables  do  not  compete  with  each  other  during  the  optimization.  The 
methodology  presented  in  the  paper  could  serve  as  a  tool  to  guide  the  design  process  during 
future  MAV  projects. 

3.  We  verified  the  inclusion  of  compliant  mechanisms  in  the  design  of  the  flapping  mechanism 
as  a  means  of  reducing  the  load  on  the  actuator.  The  mechanism  for  the  Small  Bird  MAV 
developed  by  the  University  of  Maryland  was  the  basis  of  this  work.  The  motor  torque  was 
derived  in  terms  of  rigid  body  mechanics,  and  the  design  parameters  of  the  compliant  members 
were  optimized  to  reduce  the  peak  motor  torque  during  a  flap  cycle.  Under  the  assumption 
of  the  constant  aerodynamic  load  on  the  wings,  the  results  indicated  an  89.25%  and  97.1% 
reduction  in  torque  for  the  extension  spring  and  compliant  joint  mechanisms  respectively. 

4.  We  have  also  worked  on  the  design  of  compliant  mechanisms  for  landing  gear  of  MAVs.  We 
proposed  a  jumping  mechanism  concept  that  utilizes  a  compliant  mechanism,  allowing  the 
kinematics  and  energy  storage  both  to  be  included  in  the  elastic  elements.  The  mechanism 
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weight  is  9.26g,  and  can  potentially  achieve  a  jump  height  of  lm  with  a  19.74g  payload. 
Pseudo-rigid-body  model  has  been  created  and  validated  by  Adams  simulation.  Some  further 
work  could  include  further  optimization  of  the  parameters  for  further  weight  reduction,  as  a 
factor  of  safety  of  4  is  relatively  high  and  compromise  of  the  factor  of  safety  could  be  made  for 
further  weight  savings.  Also,  different  pre-selected  dimensions  of  the  device  could  be  explored 
for  further  weight  saving. 

5.  We  have  developed  a  parameter  optimization  framework  for  determining  the  pseudo-rigid- 
body  (PRB)  model  of  cantilever  beams.  PRB  models  are  commonly  used  in  design  and 
analysis  of  compliant  mechanisms  since  they  significantly  reduce  the  number  of  degrees  of 
freedom  compared  with  the  finite  element  approach.  Also  we  proposed  a  novel  concept  of 
PRB  matrix  for  representing  topologies  of  all  PRB  models  in  a  uniform  way. 

6.  We  developed  a  theoretical  model  for  predicting  the  input  moment  required  to  actuating 
a  bistable  buckled  beam  moving  from  one  stable  to  the  other  stable  positions.  A  novel 
experimental  test  setup  was  created  for  experimental  verification  of  the  model.  The  results 
show  that  the  theoretical  model  is  able  to  predict  the  maximum  necessary  input  moment 
within  an  error  of  2.53%.  This  theoretical  model  provides  a  guideline  to  design  bistable 
compliant  mechanisms  and  actuators.  It  is  also  a  computational  tool  to  size  the  dimensions 
of  buckled  beams  for  actuating  a  specific  mechanism. 

7.  We  studied  the  effects  of  deformation  on  length  and  cross  section  of  the  beam,  and  try  to 
understand  these  phenomena  with  respect  to  some  dimensionless  parameters.  These  effects 
are  more  pronounced  for  short  beams  of  soft  materials.  A  few  PRB  models  listed  in  literature 
are  compared  against  FEA  results,  and  optimization  methods  are  used  to  determine  more 
accurate  PRB  models  for  varying  beam  geometry.  A  study  of  these  results  will  guide  the 
discussion  on  how  the  accuracy  of  the  PRB  models  changes  with  the  beam  parameters.  An 
example  of  a  compliant  mechanism  which  experiences  significant  axial  loads  will  be  used  to 
prove  the  validity  of  these  results. 

8.  We  have  developed  unified  kinetostatic  analysis  framework  for  planar  compliant  and  rigid 
body  mechanisms.  The  framework  was  written  in  MATLAB  and  is  capable  of  kinematic  and 
static  analysis  of  planar  mechanisms  with  compliant  joints  or  links.  The  simulation  results 
were  compared  with  the  Adams  software  to  test  the  validity  of  the  framework. 
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1  Classification  and  Type  Synthesis  of  Flapping  Mechanisms  for  MAVs 

Objective:  In  this  task,  we  developed  an  adjacency  matrix  based  method  for  classification  of 
transmission  mechanisms  for  flapping  wing  micro  air  vehicles.  This  method  can  be  used  to  type 
synthesis  of  flapping  mechanisms  for  MAVs. 

Methodology:  The  workspace  of  a  mechanism  is  defined  as  the  set  of  positions  that  a  workpiece 
can  reach.  The  workpiece  of  a  MAV  is  the  wing,  and  the  workspace  is  the  generated  wing  trajectory. 
The  trajectory  of  the  wing  is  responsible  for  generating  sufficient  lift  and  thrust  forces  for  flight.  The 
flapping  drive  mechanisms  are  classified  as  planar,  spherical,  or  spatial.  The  components  of  a  planar 
mechanism  all  move  in  parallel  planes,  and  utilize  lower-pair  joints.  Theoretically,  the  planar  joint 
axes  intersect  at  infinity.  A  spherical  mechanism  generates  spherically  concentric  motion  about  the 
spherical  center  of  the  mechanism.  All  the  joint  axes  must  intersect  at  the  spherical  center  point. 
Any  mechanism  which  cannot  be  classified  as  planar  or  spherical  is  a  spatial  mechanism.  However, 
planar  and  spherical  mechanisms  may  be  considered  special  cases  of  spatial  mechanisms,  due  to 
unique  geometry  and  the  orientations  of  their  joint  axes. 

Tsai  [1]  has  systematically  classified  planar  and  spatial  mechanisms  by  using  a  graph  theory. 
Each  mechanism  may  be  represented  by  a  graph  expressed  in  a  so  called  “adjacency  matrix.” 
Murphy  [2]  has  extended  this  type  synthesis  to  compliant  mechanisms  by  using  a  similar  approach. 
In  contrast  to  a  rigid  body  mechanism,  each  link  of  a  compliant  mechanism  may  be  classified  as 
rigid,  compliant,  or  ground  link.  The  connection  types  in  a  compliant  mechanism  are:  kinematic 
pair,  flexural  pivot,  or  clamped  joint.  A  compliant  mechanism  may  be  represented  by  the  so  called 
“compliant  mechanism  matrix” . 

The  adjacency  matrix  is  a  (n  x  n)  square  matrix,  where  n  de¬ 
notes  the  number  of  links  or  segments  in  the  mechanism.  For  a 
given  mechanism,  the  diagonal  elements  of  the  matrix  denote  the 
segment  type,  and  the  off-diagonal  elements  denote  the  connection 
types  between  the  segments.  The  segment  types  of  ground,  rigid 
body,  and  compliant  members  are  represented  with  a  -1,  0,  and  1, 
respectively.  The  connection  types  are  as  follows:  no  connection  is 
0,  lower  kinematic  pair  is  1,  flexural  pivot  is  2,  clamped  connection 
is  3,  gear  pair  is  4,  pulley  pair  is  5,  universal  joint  is  6,  and  scotch 
yoke  is  7.  The  composition  of  the  adjacency  matrix  is  summarized 
in  Figure  1.  The  main  focus  of  the  adjacency  matrix  is  to  cap¬ 
ture  the  driving  mechanism  of  the  MAV.  To  allow  for  universal 
comparisons  between  MAV  transmissions,  only  the  simplest  form 
of  the  mechanism  which  enables  flapping  motion  will  be  analyzed 
using  type  synthesis.  For  example  Berkeley’s  MFI  has  four  equiv¬ 
alent  linkage  mechanisms  to  ensure  stability,  type  synthesis  of  one 
of  these  mechanisms  is  sufficient.  If  the  micro  air  vehicle  had  a  relatively  simple  passive  wing 
rotation  joint  it  was  noted  and  not  included  in  the  adjacency  matrix.  However,  if  the  wing  rotation 
mechanism  exceed  a  lower  kinematic  pair,  then  a  separate  adjacency  matrix  was  developed,  i.e. 
Berkeley’s  Wing  Differential  I. 

The  most  popular  topology  in  this  survey  is  the  four-bar  linkage,  which  is  used  in  seven  of  the 
MAV  projects.  It  is  a  relatively  straightforward  design,  which  can  transform  either  a  rotational 
or  translational  input  into  a  reciprocating  flapping  motion.  The  five-bar  linkage  mechanism  is 
the  second  most  popular  with  four  of  the  MAV  projects.  Different  control  strategies  have  been 
implemented  to  utilize  this  extra  DOF,  which  will  be  discussed  in  the  next  section.  Two  six-bar 
linkages  have  been  used  in  the  selected  MAV  projects.  Both  designs  couple  two  different  planar 


Connection  Type 

None 

0 

Lower 

Kinematic 

Pair 

1 

Flexural 

Pivot 

2 

Clamped 

3 

Gear  Pair 

4 

Pulley  Pair 

5 

Universal 

Joint 

6 

Scotch 

Yoke 

7 

Segment  Type 

Ground 

-1 

Rigid  Body 

0 

Compliant 

1 

Gear 

G 

Cable 

C 

Yoke 

Y 

Figure  1:  Definition  of  segment 
types  and  connection  types 


3 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


mechanisms  to  produce  the  desired  flapping  motion.  In  the  Jumbo  Bird  MAV,  the  coupling  of  a 
crank-slider  and  a  dyad  allows  symmetry  and  stability  in  the  flapping  motion.  In  contrast,  the  MFI 
transmission  allows  amplification  of  the  linear  actuator  through  the  use  of  a  slider-crank  coupled  to 
a  four-bar  linkage.  The  last  linkage  based  system,  the  Lissajous  MAV,  is  a  double  spherical  scotch 
yoke  mechanism,  and  exhibits  a  very  cumbersome  assembly  process. 

The  Nano  Hummingbird  demonstrates  a  unique  transmission,  which  offers  advantages  and  dis¬ 
advantages  over  the  traditional  linkage-based  mechanisms.  The  Nano  Hummingbird  originally 
used  serially  coupled  four-bar  linkages,  which  transitioned  from  double-rocker  to  a  crank-rocker. 
Through  the  testing  of  prototypes,  the  design  exhibited  mechanical  wear,  radial  play  in  the  bush¬ 
ings,  decreased  efficiency,  and  failure.  The  decision  was  made  to  institute  a  double-pulley  drive,  as 
it  reduced  the  number  of  bearing  surfaces,  the  amount  of  oscillating  mass,  and  the  overall  weight 
of  the  transmission.  Another  possible  solution  may  have  been  the  implementation  of  compliant 
joints,  which  would  also  reduce  the  number  of  bearing  surfaces  and  improve  the  efficiency. 

Since  the  adjacency  matrix  was  originally  designated  for  linkage  topologies,  special  consider¬ 
ations  were  taken  for  the  double-pulley  drive.  The  driving  cable  has  a  variable  length  and  was 
treated  as  a  prismatic  joint.  The  second  cable  does  not  change  length;  therefore,  the  pulley  pairs  or 
disks  were  treated  as  if  they  had  directly  meshed  together  as  gears.  Lastly,  the  connection  between 
the  pulley  disk  and  the  variable  length  cable  was  approximated  as  a  lower  kinematic  pair. 

The  adjacency  matrices  and  the  associated  characteristic  polynomial  coefficients  are  provided  in 
Table  2  and  3.  Two  groups  of  isomorphism  were  detected  when  searching  for  equivalent  character¬ 
istic  polynomials.  The  first  isomorphism  was  identical  four-bar  linkages:  DelFly  I,  DelFly  II,  New 
DelFly  II,  and  the  LIPCA  MAV.  The  second  group  consisted  of  two  identical  five-bar  mechanisms: 
MFI  Wing  Differential  I  and  II.  It  is  interesting  to  note  that  the  HMF  and  FWMAV  are  compliant 
variations  of  the  isomorphic  four-bar  linkages,  as  seen  in  the  adjacency  matrices. 

Insect-like  kinematics  requires  control  over  three  DOF:  stroke,  pitching,  and  stroke  plane  devi¬ 
ation.  Otherwise,  a  MAV  will  approximate  insect-like  flight  through  the  use  of  additional  control 
surfaces.  Another  method  is  to  incorporate  an  extra  joint  at  the  root  of  the  wing,  which  allows  the 
wing  to  rotate  passively  along  its  long  axis.  Control  of  the  pitching  action  of  wing  has  been  noted 
in  Tables  2  and  3.  The  only  design  to  actively  control  all  three  DOF  is  the  Lissajous  MAV.  Nine 
out  of  the  fifteen  MAVs  have  a  single  DOF  transmission,  and  have  auxiliary  mechanisms  to  assist 
in  the  kinematics.  Five  transmissions  have  two  DOF,  and  utilize  different  control  strategies. 

The  first  strategy  is  to  actively  control  the  DOF,  to  generate  the  desired  wing  trajectory.  The 
MFI  is  an  example  of  this  method,  where  two  parallel  six-bar  linkage  transmissions  are  coupled  to 
a  single  wing  via  the  wing  differential.  The  parallel  six-bars  are  separately  actuated  to  control  the 
two  DOF  present  in  the  wing  differential.  The  opposite  strategy  of  passive  control  is  implemented 
in  the  PARITy  MAV.  In  the  design  iteration  from  the  HMF  to  the  PARITy,  an  extra  revolute 
joint  is  added  between  the  actuator  and  the  original  transmission.  The  idea  is  to  allow  the  passive 
balancing  of  aerodynamic  loading  on  the  wings,  which  is  similar  to  the  concept  of  an  automobile’s 
differential  during  a  turning  maneuver. 

Separate  from  the  classification  of  MAV  transmissions  is  the  classification  of  actuators,  and 
the  development  of  a  selection  database  [3].  Several  families  and  sub-types  have  been  identified 
and  categorized  by  a  number  of  characteristics  which  include:  maximum  output  strain,  maximum 
frequency,  maximum  energy  density,  efficiency,  mass,  volume,  and  many  more.  Some  of  the  major 
classes  of  actuators  are  electromagnetic,  electromechanical,  fluidic,  piezoelectric,  smart  materials, 
and  hybrid.  The  selection  of  actuator,  power  source,  control  system,  and  transmission  must  be 
integrated,  as  there  is  an  inherent  relationship  between  them. 

As  part  of  the  scalability  of  the  MAV,  the  selection  of  actuator  will  affect  the  feasibility  of  the 
design.  Scaling  actuators  from  macro  to  micro  will  result  in  varying  performance  characteristics. 
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In  some  cases,  where  a  marco  actuator  type  is  readily  available,  there  are  challenges  present  in  the 
manufacture  and  efficiency  of  the  micro  version.  In  regards  to  the  scaling  up  of  actuators  from 
micro  to  macro,  the  underlying  physics  associated  with  the  functioning  of  the  actuator  may  not 
be  possible  [4].  In  this  paper  we  will  limit  ourselves  to  the  general  categorization  of  the  type  of 
actuator  as  rotary  and  linear,  and  how  this  impacts  the  mechanism  design.  The  details  of  actuator 
selection  and  performance  are  left  to  detailed  resources. 

In  terms  of  transmission  design,  the  rotary  and  linear  actuators  require  fundamentally  different 
intermediary  mechanisms,  to  transform  the  given  input  motion  into  a  reciprocating  flapping  mo¬ 
tion.  There  are  advantages  and  disadvantages  associated  with  each  type  of  mechanism.  A  major 
difference  may  be  how  adjustable  the  kinematics  of  mechanism  is,  as  discussed  in  [5].  Adjustable 
kinematics  is  more  attainable  in  mechanism  with  linear  actuators.  Additional  DOF  can  be  added 
to  the  transmission  to  allow  the  adjustment  from  one  flight  pattern  to  another.  Rotary  based 
mechanisms  generally  have  a  constrained  output,  which  limits  the  wing  trajectory  to  a  single,  fixed 
path. 

The  majority  of  the  MAVs  included  in  this  survey  require  rotary  input.  These  nine  MAVs 
use  rotary  DC  motors,  and  the  remaining  designs  implement  linear  actuators.  The  linear  actuator 
based  transmissions  are  all  planar,  and  are  specifically  powered  by  piezoelectric  actuators.  The 
advantages  and  disadvantages  of  each  actuator  are  discussed  briefly  [6-8]. 

The  DC  motor  is  a  popular  choice  for  MAVs,  as  it  is  reliable,  versatile,  exhibits  high  efficiency, 
and  is  readily  available  in  the  market  place.  Gearboxes  are  required  to  reduce  the  motor  speed  to 
an  appropriate  level,  which  allows  the  desired  flapping  frequency  to  be  attained.  The  disadvantages 
of  DC  motors  include  the  additional  weight  of  the  gearbox,  the  additional  weight  of  a  controller  for 
brushless  motors,  and  their  limited  scalability. 

Lift  generation  depends  upon  the  wings  being  able  to  cover  a  large  range  of  flapping  angles. 
While  this  may  be  trivial  for  mechanisms  with  continuous  rotary  motion,  linear  actuators  produce 
little  strain,  and  thus  require  mechanical  amplification.  In  this  survey,  all  the  piezoelectric  (PZT) 
actuators  are  implemented  with  cantilever  boundary  conditions  and  amplified  with  some  type  of 
four-bar  linkage.  The  main  appeal  of  PZTs  is  that  they  exhibit  high  power  density  and  high 
efficiency.  A  major  disadvantage  of  PZTs  is  that  they  require  a  high  activation  voltage,  which  is 
currently  too  large  and  heavy  to  be  carried  on  board  MAVs. 

Summary  of  Main  Results  We  have  classfied  fifteen  mechanisms  in  the  literature:  ten  planar, 
two  spherical,  and  three  spatial  (Table  1).  In  general,  planar  mechanisms  have  the  advantage  of 
being  easier  to  design,  manufacture,  and  assemble.  The  major  disadvantage  of  planar  mechanisms 
for  MAV  is  the  ability  to  generate  insect-like  kinematics.  Additional  control  surfaces,  auxiliary 
mechanisms,  and  actuators  will  be  required  to  produce  a  better  approximation  of  insect-like  flight. 
Examples  of  this  are  the  DelFly  I,  New  DelFly  II,  Small  Bird,  and  Jumbo  Bird  MAVs,  which 
feature  planar  transmissions  and  tail  surfaces  [9].  Spherical  mechanisms  (i.e.  Lissajous  MAV)  are 
able  to  generate  the  figure-of-eight,  banana-shape  trajectory,  which  is  shown  in  Figure  1.  With  this 
benefit  may  come  the  complication  of  manufacture  and  assembly,  as  all  joint  axes  must  intersect  at 
the  spherical  center  of  the  mechanism  to  avoid  binding.  While  more  realistic  insect-flight  patterns 
are  possible,  additional  DOF  translates  directly  into  additional  actuators  and  more  weight.  Spatial 
mechanisms  like  spherical  mechanisms  are  able  to  generate  the  insect-like  kinematics,  but  may  offer 
more  flexibility  in  terms  of  adjustable  kinematic  parameters  (i.e.  not  limited  to  concentric  spherical 
operation).  An  example  of  this  is  the  VTOL  MAV,  which  is  able  to  generate  four  different  flight 
patterns.  When  compared  to  the  simplistic  nature  of  planar  mechanisms,  the  manufacture  and 
assembly  of  spatial  mechanisms  is  more  challenging. 

There  are  three  types  of  compliant  mechanisms,  which  will  be  considered  during  this  classi- 
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fication:  compliant  or  flexural  joints,  short  compliant  segments,  and  long  compliant  segments. 
Flexural  joints  localize  the  compliance  in  the  mechanism,  and  typically  allow  small  deflections. 
Flexural  members  distribute  compliance  in  a  mechanism  to  enhance  the  desired  transformation  of 
motion  and  forces.  Short  compliant  segments  exhibit  relatively  small  deflection,  when  compared  to 
the  large  and  possibly  non-linear  deflections  of  long  compliant  segments. 

As  part  of  this  classification,  the  transmissions  were  sorted  by  segment  type  as  follows:  seven 
rigid  body  mechanisms,  seven  compliant  joint  mechanisms,  one  short  compliant  member  mech¬ 
anism,  and  no  long  compliant  member  mechanisms.  It  is  interesting  to  note  that  only  planar 
mechanisms  feature  compliant  components;  however,  the  MFI  wing  differentials  do  feature  compli¬ 
ant  joints.  Of  the  compliant  mechanisms,  the  use  of  the  flexural  joints  dominated,  which  may  be  due 
to  their  relatively  easy  analysis  and  fabrication.  In  Section  2.2,  the  advantages  and  disadvantages 
were  discussed  in  detail. 

The  University  of  Maryland’s  Small  Bird  MAV  demonstrates  the  ideal  application  of  compliant 
mechanisms.  The  Small  Bird  uses  flexural  members  to  reduce  the  peak  motor  torque.  The  motor 
experiences  peak  torque  at  the  beginning  of  the  up-stroke  and  the  down-stroke,  due  to  the  aerody¬ 
namic  loading  of  the  wings.  When  the  crank  and  coupler  reach  over-center,  the  compliant  frame 
reaches  its  maximum  deflection.  Once  the  crank  pulls  away  from  the  maximum  point,  the  compli¬ 
ant  members  release  their  stored  energy,  and  assist  the  motor  in  overcoming  the  peak  aerodynamic 
loading  in  the  stroke  cycle.  The  compliant  frame  was  manufactured  as  a  single  piece,  using  injec¬ 
tion  molding.  This  reduces  the  use  of  traditional  joints,  decreases  the  weight  of  the  transmission, 
decreases  the  part  count,  and  improves  the  overall  efficiency  of  the  transmission.  For  a  larger  size 
MAVs,  such  as  the  Jumbo  Bird,  flexural  joints  were  used  to  avoid  long  compliant  segments,  which 
may  behave  non-linearly  and  negatively  affect  the  stability  of  the  MAV. 

The  importance  of  this  classification  lies  in  that  it  lays  out  the  foundation  of  type  synthesis  of 
flapping  wing  mechanisms.  Type  synthesis  is  the  process  of  defining  the  type  a  mechanism,  which 
is  best  suited  for  a  given  problem.  The  input  to  the  type  synthesis  process  include  (1)  type  of 
input  motion  (translational  or  rotational),  (2)  flapping  frequency  of  the  output  beating  motion,  (3) 
desired  mobility  (number  of  actuators)  and  (4)  choice  of  planar  or  spherical  or  spatial  solutions. 
This  is  typically  a  creative  stage  in  a  design  process. 

2  Optimization  of  Mechanism  Design  for  a  Flapping-Wing  Micro- Air- Vehicles 

Objective  :  Development  of  design  guidelines  regarding  the  effect  of  different  parameters  of 
flapping  wing  MAVs  on  lift  production  and  power  consumption.  Optimization  of  these  variables 
using  modeling  and  simulation  will  improve  the  performance  based  on  the  project  goals. 

Methodology:  Development  of  flapping  wing  MAVs  currently  follows  an  ad  hoc  approach,  where 
successive  designs  are  tested  until  one  is  obtained  that  satisfies  performance  requirements.  For 
example,  during  the  development  of  the  Nanohummingbird  [26],  the  team  at  AeroVironment  tried 
more  than  50  different  wing  shapes  and  structures.  A  quantitative  approach  may  be  developed 
which  can  be  undergo  an  optimization  process  so  as  to  yield  better  performance  in  the  form  of 
increased  lift,  reduced  power  consumption  or  other  key  objectives. 

Three  key  performance  characteristics  for  an  MAV  are  power  consumption,  maximum  payload 
the  vehicle  can  carry,  and  the  lift  to  power  ratio  (sometimes  replaced  by  the  thrust  to  power  ratio). 
The  important  parameters  that  affect  these  performance  characteristics  are  the  motor  characteris¬ 
tics,  the  shape  and  size  of  the  wing  and  the  mechanism  kinematics.  For  the  mechanism  performance, 
the  important  variables  are  the  dimensions  of  the  linkages.  The  mechanism  being  considered  also 
has  compliant  elements,  which  leads  to  the  definition  of  spring  constants  and  equlibrium  positions. 
By  optimizing  these  parameters,  we  can  expect  an  improvement  in  the  performance  of  the  MAV. 


6 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Table  1:  Classification  of  Flapping  Mechanisms 


Research 

Group 

Project 

Workspace 

Type 

Topology 

Rigid  Body 
or  Compliant 

Actuator 

Type 

Transmission 

Mobility 

Wing 

Rotation 

Year 

Ref. 

AeroVironment 

Nano 

Hunmmingbird 

Planar 

Double 
Pulley  Drive 

Rigid  Body 

DC  Motor 

1 

Active 

2012 

[10] 

Berkeley 

MFI 

Transmission 

Planar 

6-bar 

Compliant 

Joints 

PZT 

BiMorph 

1 

N/A 

2007 

[11] 

Wing 

Differential  I 

Spatial 

5-bar 

Compliant 

Joints 

N/A 

2 

Active 

2002 

[12] 

Wing 

Differential  II 

Spherical 

5-bar 

Compliant 

Joints 

N/A 

2 

Active 

2003 

[13] 

Warsaw/ Cranfield 
University 

Lissajous  MAV 

Spherical 

Double 
Scotch  Yoke 

Rigid  Body 

DC  Motor 

3 

Active 

2005 

[14] 

PncRa) 

VTOL  MAV 

Spatial 

4- bar 

Rigid  Body 

DC  Motor 

1 

Active 

2009 

[15] 

Delft  University 

DelFly  I 

Planar 

4- bar 

Rigid  Body 

DC  Motor 

1 

Passive 

2009 

[16] 

DelFly  II  & 
DelFly  Micro 

Spatial 

4- bar 

Rigid  Body 

DC  Motor 

2 

Passive 

2009 

[17] 

New  DelFly  II 

Planar 

4- bar 

Rigid  Body 

DC  Motor 

1 

Passive 

2010 

[18] 

Harvard 

HMF 

Planar 

4- bar 

Compliant 

Joints 

PZT 

BiMorph 

1 

Passive 

2008 

[19] 

PARITy 

Planar 

4- bar 

Compliant 

Joints 

PZT 

BiMorph 

2 

Passive 

2010 

[20,21] 

Konkuk 

University 

LIPCA  MAV 

Planar 

4- bar 

Rigid  Body 

LIPCA 

1 

Passive 

2008 

[22] 

University  of 
Delaware 

FWMAV 

Planar 

4- bar 

Compliant 

Joints 

DC  Motor 

1 

Passive 

2005 

[23] 

University  of 
Maryland 

Small  Bird 

Planar 

5-bar 

Short 

Compliant 

Segments 

DC  Motor 

2 

None 

2010 

[24] 

Jumbo  Bird 

Planar 

6-bar 

Compliant 

Joints 

DC  Motor 

1 

None 

2009 

[25] 

Compliant  mechanisms  are  used  in  the  design  since  they  reduce  problems  with  assembly  and  elim¬ 
inate  wear  and  tear  due  to  relative  motion  and  the  need  for  lubrication  [27],  and  can  also  reduce 
power  consumption.  The  base  mechanism  chosen  for  the  design  process  is  the  Jumbobird  [25] 
developed  by  University  of  Maryland. 

The  aerodynamic  forces  on  the  wing  were  calculated  using  equations  developed  by  Theodoresen 
for  the  forces  and  moments  that  act  on  an  airfoil  with  an  aileron  which  is  undergoing  flapping 
motion  in  a  stream  of  fluid  of  constant  velocity  [28].  The  lift  equation  for  the  wing  is  given  by 

L  =  pb2(7vh  +  vn ra  —  irbaa)  +  2npvbC(k)Q  (1) 

where  p  is  the  density  of  air,  b  is  half  the  chord  length  and  v  is  the  forward  velocity  of  the  wing. 
C(k)  and  Q  are  defined  in  Ref.  [28]. 

The  main  focus  of  this  work  was  to  develop  a  dynamic  model  of  the  flapping  mechanism  of 
an  MAV  so  that  the  effects  of  different  parameters  on  the  performance  of  the  MAV  can  be  easily 
studied.  The  mechanism  is  shown  in  Fig.  3.  A  modular  approach  was  used  in  the  development  of 


r 


dr 

N - H 

b 

Figure  2:  Section  of  spanwise  blade  element  of  wing  and  planform  of  wing 
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Figure  4:  Variation  of  lift  and  power  with  kinematic  parameters 


the  model  so  that  improvements  could  be  made  to  the  aerodynamic  force,  the  motor  model,  the 
compliant  mechanism  theory  and  the  optimization  without  affecting  other  parts  of  the  code.  The 
masses  and  moments  of  inertia  of  the  rigid  links  were  estimated  using  the  mechanism  dimensions. 
The  values  of  the  spring  constants  were  calculated  using  the  pseudo  rigid  body  model. 

The  MAV,  being  a  complex  system, 
has  a  large  number  of  parameters  that 
can  be  varied,  which  means  that  the 
design  space  would  extend  into  a  large 
number  of  dimension.  Since  our  goal  is 
the  optimum  design  of  the  MAV  mecha¬ 
nism,  we  looked  at  the  how  different  sys¬ 
tem  parameters  affect  the  performance 
of  the  system.  With  this  data,  we  could 
try  to  identify  the  crucial  variables  and 
work  on  optimizing  only  these  param¬ 
eters  which  would  significantly  reduce 
the  load  on  the  optimization  algorithm. 

The  design  variables  were  divided  into 
two:  (1)  kinematic  parameters  and  (2) 
system  parameters. 

The  kinematic  parameters  mainly 
refer  to  the  lengths  of  the  different  links, 
and  these  affect  the  kinematics  of  the  mechanism.  Fig.  4  shows  the  variation  of  lift  and  power  when 
each  of  the  link  lengths  was  varied  while  keeping  the  other  values  constant.  It  was  observed  that 
links  1  and  4  affected  lift  production  and  power  consumption  significantly,  whereas  link  2  affected 
only  power  consumption. 

The  system  parameters  refer  to  other  variables  in  the  system,  such  as  wing  size,  motor  char¬ 
acteristics  and  joint  stiffness.  The  important  system  parameters  which  were  incorporated  in  the 
model  were  identified  as  angle  of  attack  of  the  wing,  the  average  flapping  frequency,  the  forward 
velocity  of  the  MAV,  and  the  wing  size  which  was  determined  by  mean  chord  length  and  wingspan. 
Fig.  5  is  a  graph  showing  mean  lift  produced  with  variation  in  different  parameters  while  keeping 
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Figure  6:  Variation  of  lift  and  power  with  equilibrium  position  of  springs 


the  other  values  a  constant.  Since  the  mechanism  has  a  few  compliant  joints,  there  exists  the 
possibility  of  using  the  parameters  of  the  compliant  elements  to  achieve  the  required  system  perfor¬ 
mance.  As  explained  earlier,  the  compliant  joints  were  modeled  as  torsion  springs,  and  therefore,  it 
would  be  possible  to  vary  the  equilibrium  positions  of  the  springs.  In  reality,  this  can  be  achieved 
by  changes  in  the  fabrication  process,  such  that  the  spring  compliant  elements  are  in  a  prestressed 
state  when  assembled.  Fig.  6  describes  the  effect  of  varying  these  equilibrium  angles  on  lift  and 
power.  The  variation  in  lift  is  not  very  high,  whereas  power  varies  significantly.  It  is  clearly  bene¬ 
ficial  to  assemble  the  mechanism  such  that  certain  springs  are  always  in  a  state  of  compression  or 
expansion. 


A  series  of  optimization  runs  was  con¬ 
ducted  to  achieve  the  best  design  where  the 
variables  which  had  the  largest  influence  on 
performance  were  optimized.  A  range  of 
values  was  provided  for  each  of  the  variables 
keeping  in  mind  the  size  of  the  system  and 
other  practical  constraints.  Due  to  factors 
pertaining  to  the  accuracy  and  robustness 
of  the  aerodynamic  model,  care  was  taken 
to  ensure  that  the  flight  regime  of  the  MAV 
did  not  change  drastically.  The  aim  of  the 
optimization  was  either  l)to  produce  max¬ 
imum  lift  without  consuming  more  power, 
or  2)  to  minimize  the  power  consumption 
without  reducing  the  lift  produced. 

The  variables  used  in  the  optimization 
were  dimensions  of  two  of  the  mechanism 
links,  l\  and  Z4,  the  angle  of  attack  of  the 
wing,  cq  the  width  of  the  compliant  joints, 
u;,  and  the  initial  configuration  of  the  joints, 
defined  by  A Oq-  The  results  of  the  opti¬ 
mization  are  shown  in  Table  4. 

Fig.  7  shows  the  change  in  lift  produced  from  optimization  for  maximum  lift.  In  this  case,  the 


Figure  5:  Variation  of  lift  with  system  parameters 
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Figure  7:  Comparison  of  lift  produced  during 
flapping  motion  before  and  after  optimization 
for  maximizing  lift 


Time(s) 


Figure  8:  Comparison  of  power  consumed  dur¬ 
ing  flapping  motion  before  and  after  optimiza¬ 
tion  for  minimizing  power 


mean  power  consumption  is  not  allowed  to  exceed  the  base  value  of  0.745W.  The  graph  clearly 
demonstrates  the  increase  in  lift  with  higher  and  wider  red  peaks  compared  to  blue  peaks.  Fig.  8 
shows  the  possible  reduction  in  power  consumption  when  35g  of  lift  is  required.  Since  the  negative 
power  values  are  ignored,  it  is  seen  that  the  positive  peaks  are  smaller  than  before  optimization, 
and  should  result  in  lower  power  consumption. 

Close  inspection  of  Table  4  suggests  that  increasing  the  constraint  on  power  consumption  does 
not  provide  a  substantial  increase  in  lift  produced.  For  example,  changing  the  maximum  power 
consumption  from  0.745W  to  0.9  W  only  increases  lift  from  48.95g  to  54.72g.  However,  reducing 
the  mass  of  the  system,  that  is,  the  minimum  lift  required  from  35.12g  to  30g  reduces  the  power 
consumption  from  0.4044W  to  0.2151W.  This  means  that  it  is  more  beneficial  to  focus  on  reducing 
the  mass  of  the  system  compared  to  using  a  more  powerful  motor.  It  was  also  noticed  during  the 
optimization  that  it  is  beneficial  to  have  the  springs  in  the  system  in  such  a  configuration  that  they 
aid  the  downward  motion  of  the  wing  (during  which  the  lift  is  produced).  The  optimization  moves 
towards  the  region  where  the  equilibrium  position  of  the  springs  is  at  the  bottom  of  the  flapping 
motion. 

Summary  of  Main  Results  :  The  dynamic  model  is  a  useful  tool  to  study  the  performance 
of  the  system  with  respect  to  different  design  parameters.  The  sensitivity  study  provides  valuable 
insight  into  contribution  of  individual  parameters  towards  the  generation  of  lift  as  well  as  power 
consumption.  Once  these  parameters  are  identified,  the  design  can  be  optimized  for  required 
performance  goals. 

A  salient  feature  of  a  majority  of  the  results  was  that  the  optimal  values  were  usually  at  the 
upper  or  lower  bounds  of  the  variables,  which  means  that  improvements  in  the  performance  of 
the  MAV  appear  to  be  limited  only  by  fabrication  and  size  constraints.  However,  some  of  the 
assumptions  involved  in  the  development  of  the  model  must  be  kept  in  mind  while  attempting 
to  change  the  parameters.  One  of  the  significant  limitations  is  the  reliability  of  the  aerodynamic 
model  in  accurately  predicting  the  forces  on  the  wing. 

Depending  on  the  compliant  joint,  the  pseudo-rigid-body  model  can  be  altered  for  more  accurate 
results.  Some  accuracy  issues  can  be  overcome  by  developing  a  better  system  model,  which  may 
require  higher  computing  ability.  The  aerodynamic  model  itself  can  be  improved  by  incorporating 
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all  the  dominant  aerodynamic  effects.  Considering  how  complex  such  a  force  model  would  be  if 
derived  from  basic  principles,  a  heuristic  or  surrogate  model  would  be  much  more  compatible. 

A  summary  of  the  design  guidelines  is  given  below. 

1.  If  compliant  joints  are  present,  they  should  be  assembled  such  that  the  stiffness  aids  the 
downward  motion  of  the  wing.  This  can  also  be  achieved  by  using  regular  springs. 

2.  The  focus  of  the  design  process  should  be  on  reducing  the  mass  of  the  system  rather  than 
increasing  the  power  available  to  it. 

3.  Increasing  the  flap  angle  leads  to  simultaneous  increase  in  lift  production  and  power  con¬ 
sumption.  For  the  system  described  here,  the  increase  in  power  is  greater  than  the  increase 
in  lift. 

The  approach  described  provides  a  set  of  guidelines  for  design  of  a  flapping  wing  MAV  for 
suitable  dynamic  performance.  This  method  can  be  altered  to  fit  the  application  in  question,  but 
they  are  general  principles  that  can  be  used  for  any  dynamic  compliant  mechanism.  With  the  choice 
of  a  suitable  aerodynamic  model,  the  external  forces  on  the  system  can  be  determined  with  minimal 
computation.  A  simple  dynamic  model  can  be  used  to  approximate  the  performance  of  the  system, 
and  the  sensitivity  study  determines  the  important  design  parameters.  Using  this  information,  an 
optimization  algorithm  can  be  utilized  to  arrive  at  the  mechanism  design  that  leads  to  optimal 
performance. 

3  A  3-spring  Pseudo  Rigid  Body  Model  for  Soft  Joints  with  Significant  Exten¬ 
sion  Effects 

Objective:  The  goal  is  to  derive  a  rigid  body  approximation  that  can  capture  the  extension 
effects  in  similar  compliant  joints.  We  have  derived  a  3-spring  revolute-prismatic-revolute  pseudo- 
rigid-body  model  for  short  beams  used  in  soft  joints  made  of  elastomer  material.  These  elastomer 
joints  are  often  used  in  compliant  transmission  mechanisms  for  flapping  wing  MAVs. 

Methodology:  Compliant  mechanisms  are  those  in  which  conventional  rigid-body  elements  such 
as  pin  joints  are  replaced  by  deformable  members  or  flexures.  Methods  for  analysis  and  design  of 
compliant  mechanisms  have  improved  significantly  over  the  last  couple  of  decades.  This  work  deals 
with  a  technique  that  uses  the  pseudo  rigid  body  (PRB)  model,  most  of  the  groundwork  for  which 
was  laid  down  by  Howell  and  Midha  [29,30].  Since  the  behavior  of  deformable  members  is  difficult 
to  introduce  into  the  kinematic  study  of  a  mechanism,  they  are  converted  to  rigid  body  links  for 
the  ease  of  analysis.  For  instance,  a  flexure  pivot  may  be  replaced  by  a  revolute  joint  and  a  torsion 
spring. 

Many  PRB  models  have  been  developed,  each  more  accurate  than  the  previous  one,  but  with 
accuracy  comes  the  cost  of  complexity.  For  beam  type  compliant  elements,  Dado  [31]  developed  a 
variable  parametric  model  and  Su  [32]  came  up  with  a  model  with  3  revolute  joints  .  However,  most 
PRB  models  only  account  for  the  transverse  deformation  of  the  beams,  while  the  axial  deformation 
of  the  beams  is  ignored.  In  this  work,  the  model  presented  will  also  take  axial  deformation  into 
consideration.  Vogtmann  et  al.  [33]  developed  a  similar  model  for  elastomeric  joints  to  account  for 
extension  and  torsion  effects,  demonstrating  the  need  to  include  extension  springs  for  soft  joints. 
It  provided  a  quick  derivation  of  the  PRB  parameters,  and  was  shown  to  be  accurate  for  specific 
types  of  joints.  The  approach  in  this  paper  is  more  general  and  detailed,  and  introduces  a  rigid 
segment  at  either  end  of  the  beam  which  is  shown  to  affect  the  error  in  the  model. 

The  PRB  model  has  to  be  compared  with  an  analytical  or  an  FEA  or  a  physical  model  to 
determine  the  values  of  the  parameters.  For  this  case,  a  beam  model  approximation  based  on 
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Figure  10:  The  3-spring  PRB  model  after  de¬ 
formation 
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Figure  11:  Optimal  PRB  parameters  for  RPR 
model  as  functions  of  length  to  height  ratio 


Timoshenko  beam  theory  [34,  35]  has  been  developed.  The  model  is  analyzed  over  a  large  range 
of  general  loads  and  a  comparison  is  made  between  the  Timoshenko  beam  model  and  PRB  model. 
An  optimization  process  is  carried  out  to  determine  the  PRB  model  parameters  that  give  the  least 
error.  Figure  9  shows  a  few  soft  joints  using  the  SDM  technique  at  the  Design,  Innovation  and 
Simulation  Lab  (DISL).  The  flexible  parts  are  made  from  a  rubber-like  material  of  very  low  modulus 
of  elasticity.  During  testing,  it  was  seen  that  due  to  the  low  elastic  modulus  of  the  material,  the 
flexures  also  undergo  axial  deformation  and  Poisson’s  effect.  A  beam  model  was  developed  that 
takes  these  effects  into  account  and  used  it  to  derive  a  pseudo-rigid-body  model  for  computational 
purposes. 

In  order  to  account  for  the  various  be¬ 
havioral  characteristics  of  the  soft  joint,  a 
pseudo-rigid-body  model  of  three  segments 
was  developed.  The  model  is  a  revolute- 
prismatic-revolute  (RPR)  serial  chain.  A 
schematic  of  the  same  is  shown  in  Fig.  10. 

The  schematic  is  for  a  beam  of  unit  length. 

The  base  model  is  symmetric,  i.e.  the  pa¬ 
rameters  are  the  same  at  either  end.  There 
are  two  rigid  beams  of  length  7.  An  ex¬ 
tension  spring  of  stiffness  Kex  is  attached 
to  both  the  rigid  members  through  revolute 
joints.  Torsion  springs  are  present  at  both 
the  revolute  joints,  of  spring  constant  Kq. 

The  extension  spring  is  not  capable  of  bend¬ 
ing,  as  if  a  prismatic  joint  is  present  between 
the  two  revolute  joints.  The  torsion  springs 
account  for  the  bending  stiffness  of  the  beam,  and  the  extension  spring  represents  the  axial  stiffness. 
The  undeflected  length  of  the  extension  spring  is  Zo- 

A  modified  Timoshenko  beam  approximation  was  used  to  determine  the  deflection  of  the  beam 
under  the  action  of  tip  loads,  including  bending,  shear,  extension  and  Poisson’s  effect.  The  results 
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Figure  9:  Joints  fabricated  using  shape  deposition 
manufacturing  at  Design,  Innovation  and  Simulation 
Lab  at  OSU. 
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Figure  12:  Variation  of  error  in  PRB  model  with  loads  (for  the  joint  with  L/h  —  2).  ax  and  ay  are  tip 
forces  in  the  axial  and  transverse  directions.  /3  is  the  bending  moment  at  the  beam  tip. 


from  this  method  were  used  to  calculate  the  optimal  values  of  the  PRB  parameters  over  a  range  of 
loads. 

The  extension  effects  in  the  soft  joints  are  expected  to  change  with  respect  to  the  ratio  of  length 
to  width  for  the  compliant  joint.  Nine  values  of  the  length  to  height  ratio  (for  beams  of  rectangular 
cross-section)  were  considered,  with  the  L/h  ratio  varying  between  2  and  10.  The  optimal  values  of 
the  PRB  parameters  obtained  for  each  case  are  shown  in  Fig.  11,  with  the  curve  fitting  functions 
given  below. 

-0.028lg  +  0.2266^  -  0.2698 
0.0018^  +  0.1792 

ke  =  -0.0017f  +  2.0220 
h 

kex  =  -0.0942  f  +  1.3002 
h 

Figure  12  shows  the  details  on  the  variation  of  error  with  the  loads  viewed  as  contour  slice  plots. 
A  thorough  reading  of  the  plots  indicates  that  the  error  is  very  low  (below  0.02)  for  a  large  range  of 
loads,  and  the  only  regions  of  high  error  are  when  all  the  loads  work  towards  deflecting  the  beam 
in  the  same  direction,  i.e.  ax  =  —  1  and  ay  =  ft  =  ±1.  The  error  plots  clearly  indicate  that  errors 
creep  into  the  model  only  for  very  large  loads,  and  it  works  extremely  well  for  tip  deflections  as 
high  as  70°.  A  sensitivity  analysis  was  also  performed  near  the  optimal  point  to  determine  the 
significance  of  each  PRB  parameter  on  the  error.  The  results  of  the  study  for  the  shortest  beam 
(L/h  =  2)  are  presented  in  Fig.  13.  A  reading  of  the  vertical  axis  of  the  graphs  also  indicates  that 
the  change  in  error  due  to  kex  is  much  less  than  the  change  in  error  due  to  either  k@  or  7.  Thus, 
the  results  are  more  sensitive  to  changes  in  the  latter  two  parameters. 

Summary  of  Main  Results:  The  pseudo-rigid-body  model  approximation  described  in  the 
form  of  a  revolute-prismatic-revolute  (RPR)  member  is  capable  of  representing  the  bending  and 
extension  stiffness  of  the  compliant  soft  joint.  Due  to  the  combination  of  bending  and  extension 
effects,  the  PRB  parameters  are  dependent  on  the  aspect  ratio  of  the  beam,  but  still  independent 
of  the  modulus  of  elasticity  and  can  be  used  for  different  dimension  scales. 

The  PRB  model  presented  here  is  limited  by  the  effectiveness  of  the  beam  theory  used  to 


2  <  L  <  4 

4  <  \  <  10 


13 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


calculate  the  force-deflection  relationship.  With  a  better  understanding  of  the  material  behavior 
and  better  solution  techniques,  it  is  possible  to  improve  the  characteristics  of  the  model.  Currently, 
it  is  limited  to  the  linear  elastic  regime  of  the  material.  The  deformation  of  the  members  is  also 
restricted  to  bending  and  extension,  without  considering  distortion  of  the  cross-section. 

The  integration  of  extension  effects  into  the  analysis  of  compliant  members  and  the  PRB  model, 
the  optimization  approach  and  the  fitting  functions  of  the  PRB  parameters  are  the  main  contri¬ 
butions  of  this  work.  The  simplicity  of  the  PRB  models  when  compared  to  beam  models  or  FEA 
makes  them  ideal  candidates  for  design  and  analysis  of  compliant  mechanisms,  for  computational 
reasons  or  when  using  an  intuitive  design  approach.  By  reducing  the  error  and  developing  a  general 
framework,  it  is  possible  to  create  a  design  approach  for  compliant  mechanisms  using  PRB  models. 

4  A  Parameter  Optimization  Framework  for  Determining  the  Pseudo-Rigid- 
Body  Model  of  Cantilever  Beams 

Objective  :  Development  of  a  framework  for  determining  the  optimal  pseudo-rigid-body  (PRB) 
model  of  2D  cantilever  beams.  This  method  can  be  used  for  deriving  numerical  approximations  for 
compliant  elements  for  various  applications. 

Methodology:  The  design  of  compliant 
mechanisms  can  be  a  difficult  process  due 
to  the  nonlinearities  arising  out  of  their  de¬ 
flections.  This  places  a  great  emphasis  on 
the  need  to  improve  the  methods  of  anal¬ 
ysis.  With  a  better  understanding  of  be¬ 
havior  of  compliant  elements  and  enhanc¬ 
ing  the  tools  for  analysis,  there  will  be  a 
significant  advancement  in  the  use  of  these 
mechanisms.  The  approach  used  in  this 
work  replaces  flexible  elements  with  pseudo- 
rigid-body  (PRB)  models,  which  simplifies 
the  statics  and  kinematics  equations.  Since 
rigid-body  kinematics  are  well  studied,  the 
PRB  approach  is  usually  more  intuitive  for 
the  purpose  of  analysis  and  design.  How¬ 
ever,  care  needs  be  taken  to  ensure  that  the 
model  is  accurate  in  replicating  the  behav¬ 
ior  of  the  flexible  members. 

Apart  from  the  articles  mentioned  in 
Sec.  3,  multiple  research  groups  have  fo¬ 
cused  on  PRB  models  over  the  past  decade  or  so.  Chen  et  al.  [36]  used  a  particle  swarm  optimizer 
to  suggest  better  values  for  the  parameters  of  the  3R  model  derived  by  Su  [32].  But  this  model 
may  be  too  complex  for  simple  elements,  leading  to  increased  computational  cost.  Recently,  Yu  et 
al.  [37]  proposed  a  2R  model  that  reduces  the  complexity  of  the  3R  PRB  model  while  maintaining  a 
similar  level  of  accuracy.  Another  method  for  using  PRB  models  with  revolute  joints  was  suggested 
by  Pei  et  al.  [38].  Some  models  also  use  prismatic  joints  with  extension  springs,  such  as  the  one 
presented  by  Saxena  and  Kramer  [39]  or  Vogtmann  et  al  [40]. 

One  of  the  major  drawbacks  of  most  of  these  models  is  the  dependence  of  the  parameters  on 
the  loading  direction.  Although  this  is  useful  for  static  analysis,  where  the  loads  and  boundary 
conditions  are  known,  it  makes  these  models  difficult  to  use  for  dynamic  simulations  or  mechanism 


Figure  13:  Sensitivity  of  individual  PRB  parameters 
near  the  optimal  point  for  RPR  model 
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Figure  14:  General  PRB  model  with  3  elements  with  its  PRB  matrix 


synthesis.  Load  independent  models  are  more  useful  for  these  cases.  Since  compliant  mechanisms 
usually  have  a  limited  range  of  motion,  the  models  only  have  to  be  accurate  within  the  expected 
range  of  loading  or  deformation.  It  may  also  be  important  to  have  model  parameters  based  on  the 
type  of  deformation  expected. 

The  required  accuracy  and  kinematics  of  the  PRB  model  also  depend  on  the  application.  It  may 
be  useful  to  develop  a  model  on  a  case  to  case  basis  after  studying  expected  deformation,  required 
accuracy  and  other  application  requirements.  With  this  in  mind,  a  strong  argument  can  be  made 
for  a  general  PRB  model  that  can  be  adapted  for  each  case.  The  model  can  be  altered  to  be  simple 
or  complex  based  on  the  demands  of  the  application.  Since  the  design  process  generally  involves 
multiple  evaluations,  a  simple  model  with  minimum  parameters  would  be  suitable  for  this.  For 
analysis,  the  user  may  be  better  served  having  a  complex  model  that  captures  all  the  characteristic 
behavior  of  the  compliant  element. 

Pseudo-rigid-body  models  serve  as  replacements  for  compliant  elements.  This  means  that  they 
should  have  properties  similar  to  those  of  the  original  elements  they  replace.  Compliant  members 
can  undergo  deformation,  but  also  provide  a  resistance  to  deformation  due  to  the  intrinsic  material 
properties.  In  order  to  capture  these  properties,  most  PRB  models  consist  of  joints,  which  define 
the  degree  of  freedom,  and  springs,  which  represent  the  stiffness  of  the  material.  If  the  fundamental 
elements  that  can  represent  these  characteristics  are  identified,  it  is  possible  to  develop  any  PRB 
model  as  a  combination  of  these  elements.  This  is  the  basic  idea  behind  the  general  PRB  model. 

Since  most  of  the  PRB  models  currently  in  literature  only  involve  revolute  joints  with  torsion 
springs  and  prismatic  joints  with  extension  springs  [41],  the  scope  of  this  discussion  will  be  limited 
to  similar  segments.  In  theory,  it  is  possible  to  add  any  segment  that  is  representative  of  any  kind 
of  deformation  or  stiffness.  Additionally,  only  springs  with  constant  stiffness  will  be  considered, 
although  nonlinear  stiffness  functions  can  also  be  expressed  using  this  method.  Figure  14  shows  a 
schematic  diagram  of  a  general  PRB  model  with  3  elements. 

One  important  aspect  of  the  general  PRB  model  is  the  number  of  parameters  needed  to  define 
it.  In  the  basic  form  discussed  in  this  paper,  each  segment  has  3  parameters,  which  translates  to  9 
parameters  for  a  model  with  3  segments.  A  convenient  method  of  representation  is  necessary  for 
these  parameters.  To  this  end,  a  simple  tool  which  called  the  PRB  matrix  was  defined.  This  matrix 
consists  of  the  values  of  the  different  parameters,  with  the  number  of  rows  equal  to  the  number 
of  segments  in  the  model.  Each  row  contains  the  values  of  k$,  kex  and  7  for  each  segment.  The 
symbol  is  used  to  represent  the  PRB  matrix.  A  representation  of  a  3x3  PRB  matrix  is  shown  on 
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the  right  side  of  Fig.  14. 

The  estimation  of  PRB  parameters  can  be  done  using  a  variety  of  methods.  Most  of  the  early 
work  on  PRB  models  focused  on  understanding  the  behavior  of  the  beam,  and  deriving  the  PRB 
model  and  the  values  of  the  parameters  based  on  observed  results  [42,43].  While  this  works  very 
well  on  a  case  to  case  basis,  it  is  very  difficult  to  use  in  an  algorithm  to  find  the  most  optimal  PRB 
model.  A  few  new  papers  have  talked  about  using  different  mathematical  techniques  for  arriving 
at  the  optimal  PRB  parameter  values  [36].  An  effective  method  is  to  perform  a  straightforward 
optimization  on  the  parameters  with  the  objective  of  trying  to  minimize  the  error  in  the  model 
compared  to  results  from  established  methods  such  as  beam  theory,  FEA  or  experiments.  In  this 
work,  the  optimal  PRB  parameters  were  calculated  using  a  local  optimization  routine  to  minimize 
the  error  between  the  PRB  model  and  the  beam  theory  results  over  a  large  number  of  loading  cases. 

Tables  5,  6  and  7  show  various  possible  topologies  of  the  PRB  model  using  combinations  of 
the  three  basic  types  of  segments.  The  tables  also  list  the  degrees  of  freedom  of  the  model  and 
the  number  of  independent  PRB  parameters.  The  error  value  represents  the  mean  error  over  the 
range  of  loads  for  the  optimized  version  of  each  model.  The  optimal  values  of  the  parameters  are 
also  listed.  In  cases  when  similar  PRB  models  have  already  been  listed  in  literature,  the  reference 
number  has  been  provided  next  to  the  name. 

It  is  possible  to  use  a  simple  algorithm  to  choose  the  best  PRB  model  for  a  given  application. 
First,  the  deflection  of  the  relevant  compliant  member  must  be  determined  using  beam  theory, 
FEA  or  experiments  for  various  loading  cases.  An  error  calculation  function  must  be  developed 
to  determine  the  error  over  all  loading  cases  for  any  given  PRB  matrix.  Using  tables  5,  6  and  7 
as  guides,  an  optimization  algorithm  can  be  used  to  run  through  various  PRB  models  to  find  the 
most  optimal  one  for  the  application. 

Summary  of  Main  Results:  The  definition  of  the  PRB  matrix,  which  allows  the  uniform 
representation  of  many  pseudo-rigid-body  models,  is  an  important  contribution.  It  can  be  used 
to  represent  varying  topologies  of  PRB  models  and  also  the  different  elements  in  them.  The 
development  of  the  optimization  framework,  which  facilitates  the  user  to  find  the  best  PRB  model 
with  optimal  parameters  for  a  given  application,  is  the  most  important  contribution.  In  this  work, 
a  direct  optimization  method  has  been  used  to  calculate  the  optimal  values  in  the  PRB  matrix.  The 
definition  of  the  error  function  is  subjective,  and  may  be  changed  depending  on  the  application. 

The  general  PRB  model  could  serve  to  be  a  useful  tool  in  the  analysis,  and  possibly  synthesis 
of  compliant  mechanisms  especially  in  the  early  design  stages.  It  eliminates  the  need  to  determine 
the  best  PRB  model  by  an  ad  hoc  approach.  Moreover,  this  procedure  can  be  easily  extended  to 
a  variety  of  compliant  elements.  Also  the  PRB  matrix  will  be  implemented  into  a  computational 
design  software  for  compliant  mechanisms  currently  being  developed  at  The  Ohio  State  University. 

5  DAS-2D:  A  Concept  Design  Tool  for  Planar  Compliant  Mechanisms 

Objective:  Currently  design  and  analysis  of  compliant  mechanisms  rely  on  several  commercial 
dynamics  and  finite  element  simulation  tools.  However  these  tools  do  not  implement  the  most 
recently  developed  theories  in  compliant  mechanism  research.  In  this  task,  we  present  Compli¬ 
ant  DAS-2D  (Design,  Analysis  and  Synthesis),  a  conceptual  design  tool  which  integrates  the  re¬ 
cently  developed  pseudo-rigid-body  models  and  kinetostatic  analysis/synthesis  theories  for  compli¬ 
ant  mechanisms. 

Methodology:  In  terms  of  computer-aided  design  of  mechanisms,  there  has  been  a  number 
of  static  solvers  that  have  been  developed  throughout  the  years  and  KINMAC  (KINematics  of 
MAChinary)  and  STATMAC  (STATics  of  MAChinary)  by  Paul  [44]  are  maybe  the  first  examples 
of  the  kinematic  and  static  analysis  tools.  SAM  [45]  is  a  commercial  software  capable  of  kinematic 
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and  static  analysis  limited  to  rigid-body  linkages.  More  recently  CoMeT  [46]  [47]  was  developed  for 
static  analysis  of  compliant  mechanisms  which  is  not  limited  to  just  the  planer  cases.  However  direct 
stiffness  method  was  employed  in  the  software  which  is  only  accurate  for  small  deflections  of  elastic 
members  limiting  the  capabilities  of  the  software.  Geometrically  Exact  Beam  Theory  (GEBT)  [48] 
is  a  FEM  tool  for  static  analysis  of  structures  but  it  is  not  possible  to  analyze  mechanisms  that 
contain  moving  joints.  Similarly,  SPACAR  [49]  is  a  software  package  for  dynamic  modeling  and 
analysis  of  multibody  flexible  systems  with  finite  element  analysis.  WorkingModel  2D  is  a  commercial 
software  for  multi-body  dynamic  simulation  of  planar  physical  systems.  It  has  been  often  used  in 
design  and  simulation  of  planar  mechanisms.  Matlab’s  SimMechanics  toolbox  is  a  powerful  multi¬ 
body  dynamics  solver  however  the  design  process  with  block  diagrams  is  not  entirely  intuitive  and 
the  toolbox  lacks  of  proper  interactivity  between  the  user  and  the  simulation.  Yue  et  al.  [50]  have 
developed  virtual  reality  user  interfaces  that  communicate  with  SimMechanics  solver  for  interactive 
design  of  planar  mechanisms.  Adams  [51]  is  a  commercial  software  capable  of  dynamic  simulation 
of  planar  mechanisms. 

Given  numerous  advances  in  compliant  mechanism  theories  in  the  past  two  decades,  they  have 
not  been  implemented  into  any  design  software  that  is  dedicated  to  design  of  compliant  mechanisms. 
Here  we  aim  to  fill  this  gap  by  developing  a  conceptual  design  software  that  integrates  kinetostatic 
analysis  and  synthesis  theory  for  design  of  compliant  mechanisms.  Kinematics  of  rigid-body  mech¬ 
anisms  will  also  be  a  side  product  of  this  software  since  they  are  considered  as  special  subset  of 
compliant  mechanisms.  The  software  was  developed  with  the  graphic  user  interface  in  MATLAB 
to  take  advantage  of  the  built-in  functions  of  nonlinear  equation  solvers  and  optimization  routines. 

It  is  well  known  that  analysis  of  compliant  mechanisms  involves  simultaneously  solving  a  set  of 
kinematic  equations  coupled  with  static  force  equilibrium  equations,  which  are  called  kinetostatic 
(kinematic  and  static)  equations.  Complex  number  method  is  one  of  the  most  commonly  used 
method  in  kinematic  analysis: 

—  re10  —  r(cos#  +  i  sin  9)  (2) 

where  i  —  y/—l  is  the  complex  number  unit.  In  general,  the  vector  loop  equation  for  a  mechanism 
with  l  independent  closed  loops  can  be  written  in  the  complex  form: 

n 

Y^  =  0,  1  <k<l  (3) 

2=1 

where  coefficient  Cki  can  be  -1,  0  or  1. 

Graph  theory  is  often  employed  in  representing  mechanisms  [52-54].  This  representation  is 
widely  used  in  classification  and  type  synthesis  of  rigid-body  mechanisms.  The  two  different  graph¬ 
ical  representations  of  the  same  slider-crank  mechanism  are  showed  in  Fig.  15.  At  the  top,  Freuden- 
stein  and  Maki  [55]  give  an  example  of  the  common  practice  in  which  vertices  represent  links  and 
the  edges  for  kinematic  joints.  Since  rigidly  connected  vertices  are  represented  with  a  single  vertex, 
this  representation  is  not  well  suited  for  implementation.  Also  only  binary  joints  can  be  handled 
with  this  representation  without  breaking  multiple  joints  into  binary  joints.  Therefore  a  new  repre¬ 
sentation  is  adapted  over  the  traditional  method.  In  this  new  representation,  vertices  represent  the 
nodes  (or  the  joints)  and  the  edges  denote  the  links  between  the  joints.  This  new  representation  is 
parallel  with  how  individual  links  are  stored  in  the  software. 

The  main  challenge  of  kinematic  analysis  is  finding  independent  loops  for  formulating  kinematic 
constraint  equations.  In  the  graph  theory,  the  topology  of  any  mechanisms  is  represented  in  the 
adjacency  matrix.  And  the  number  of  independent  kinematic  loops  can  be  calculated  with  the 
Euler’s  formula.  Once  the  adjacency  matrix  and  the  number  of  independent  kinematic  loops  are 
determined,  graph  theory  can  be  used  to  find  the  cycle  bases  (independent  loops). 
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Figure  15:  A  slider-crank  mechanism  its  the  graph  theory  representations:  the  classical  graph  repre¬ 
sentation  (top)  and  the  new  graphical  representation  implemented  in  Compliant  DAS-2D. 

After  all  independent  kinematic  loops  are  found,  nonlinear  kinematic  equations  of  the  mecha¬ 
nism  can  be  readily  derived  as  the  following.  The  kinematic  equation  for  the  zth  link  can  be  written 
as: 

—  Zio  +  A i  x  fi  (xi),  (4) 

where  Z^  and  A i  are  constant  design  parameters  for  the  link  and  ^  is  a  function  for  which  input 
is  the  unknown  parameter(s)  needed  to  define  the  link.  As  an  example,  for  a  binary  link  with  two 
pin  joints,  the  term  is  zero,  \  —  is  length  of  the  link.  And  /  =  (cos^,sin^)T  where  Oi  is 
the  link  angle. 

Following  this  notation,  for  a  mechanism  with  n  links  and  l  independent  loops,  the  kinematic 
constraint  equations  can  be  systematically  derived  as  n  vector  equations  or  2 n  scalar  equations: 

J2Cki^i  =  {  q  | ,  1  <  k  <1  (5) 

where  coefficient  Cki  can  be  -1,  0  or  1. 

By  applying  Eq.(5),  kinematic  equations  will  be  formed  and  they  can  be  solved  simultaneously 
using  a  nonlinear  equation  solver  with  the  number  of  inputs  equal  or  less  then  the  degrees  of  freedom 
of  the  mechanism. 

Compliant  mechanisms  have  at  least  one  compliant  link  that  can  be  deformed  under  external 
forces.  PRB  models  have  been  widely  used  in  the  analysis  and  synthesis  of  compliant  mechanisms. 
In  this  approach,  n  rigid  segments  are  joined  with  torsion  and  or  extension  springs.  Fig.  16  shows 
example  PRB  model  of  a  cantilever  beam.  Several  PRB  models  were  developed  to  analyze  elastic 
members. 


Figure  16:  The  PRB-3R  model  (left)  and  the  finite  segment  model  (right)  for  compliant  beams. 

Recently,  Venkiteswaran  and  Su  [56]  proposed  a  uniform  way  of  representing  different  PRB 
models  with  a  matrix  called  a  PRB  matrix: 


n  = 


7l  ^exl  ^91 


Ifn  kexn  ^ On 


(6) 
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where  each  row  represents  a  different  segment  in  the  flexible  beam.  Each  row  consists  of  length 
of  the  segment (7),  torsion  spring (kex)  and  the  flexibility  of  the  spring(fc^ri).  This  representation  of 
compliant  beams  is  adopted  in  the  software  and  each  individual  flexible  beam  can  be  represented 
with  a  predefined  but  customizable  PRB  matrix. 

Kinetostatic  analysis  is  to  determine  the  deflection  of  a  compliant  mechanism  upon  a  given 
load  or  vice  versa.  To  derive  kinetostatic  equations  (coupled  kinematic  and  static  equilibrium), 
several  methods  for  derivation  of  kinetostatic  formulation  exist.  Direct  stiffness  approach  [57-59] 
can  be  used  to  perform  linear  static  analysis.  This  method  is  currently  being  used  in  static  analysis 
of  rigid-body  (SAM  [45])  and  compliant  mechanisms  (CoMeT  [47]). Although  it  has  advantage  of 
simple  formulation,  the  stiffness  approach  is  based  on  linear  models,  hence  not  accurate  for  large 
deflections. 

On  the  other  hand,  energy  methods  [60]  can  be  used  for  nonlinear  static  analysis.  Virtual  work 
principle  was  employed  to  calculate  deformation  of  parallel  mechanisms  [61,62].  Optimization  of 
potential  energy  has  been  used  for  a  wide  variety  of  applications  from  molecular  mechanics  [63] 
to  human  walking  [64].  In  this  task,  we  implemented  the  energy  based  approach  for  formulation 
of  kinetostatic  equations.  A  stationary  potential  energy  is  the  necessary  and  sufficient  condition 
for  equilibrium  of  a  conservative  system.  This  means  that  equilibrium  position  (stationary)  point 
can  be  obtained  from  minimization  of  the  potential  energy  function  and  work  for  external  forces 
or  moments.  Furthermore,  this  optimization  is  also  subject  to  the  kinematic  constraint  defined  in 
Eq.(3).  In  summary,  the  optimization  problem  to  be  solved  is  mathematically  formulated  as: 


min  f  — 

t i 


subject  to  g  — 


2=1 


(7) 


where  (fc^,^)  equals  to  (ke,x)  and  (k$,9)  for  linear  and  torsion  springs,  respectively. 

Summary  of  Main  Results  This  task  resulted  in  DAS  2D  [65],  a  conceptual  design  tool  for  pla¬ 
nar  compliant  mechanisms.  The  compliant  mechanism  design  software  is  implemented  in  MATLAB 
using  object  oriented  programming  by  taking  advantage  of  the  built-in  nonlinear  equation  solver, 
optimizer  and  plotting  tools.  Also  graphical  user  interfaces  have  been  implemented  to  facilitate  the 
design  process.  We  developed  five  different  analysis  modules  (Fig.  17)  capable  of  kinematic  and 
kinetostatic  analysis  of  mechanisms. 

Static  force  analysis  is  to  determine  the  relationship  between  the  deflection  of  a  compliant 
mechanism  and  the  external  load  applied  upon.  The  module  starts  with  converting  compliant  links 
into  a  series  of  rigid-body  links  using  pseudo-rigid-body  model  described  earlier.  Then,  kinematic 
equations  are  derived  based  on  the  independent  kinematic  loops.  Afterwards,  a  breadth-first  search 
algorithm  is  used  to  determine  the  deflection  of  the  point  where  the  external  force  is  applied.  Finally 
the  optimizer  in  Eq.(7)  is  employed  in  the  last  optimization  step  for  the  static  analysis. 

Opposite  to  the  static  force  analysis,  the  distance  analysis  is  to  determine  the  required  load(s) 
that  will  result  in  a  prescribed  deformation.  The  desired  deformation  can  be  translation  of  a  node 
or  rotation  of  a  link  and  the  output  of  the  module  will  be  the  magnitudes  of  the  loads  exerting  at 
a  given  position. 

Mechanical  advantage  is  defined  as  the  ratio  of  the  output  force/moment  over  the  input  force/moment. 
It  is  one  of  the  most  important  design  criterion  in  compliant  mechanisms.  Mechanical  advantage 
analysis  is  to  find  the  balancing  load  for  a  given  input  load. 
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a)Kinematic  Analysis 


b)Static  Analysis 


c)Distance  Analysis 


d)Mechanical  Advantage 


e)Energy  Analysis 


Figure  IT:  Five  different  analysis  modules  are  shown. 

The  energy/torque  vs.  the  input  load  curve  gives  the  designer  a  visualized  way  to  evaluate  the 
quality  of  a  compliant  mechanism.  The  energy  analysis  module  plots  the  distance  versus  load  and 
energy  graph  until  a  desired  deflection. 

6  DAS-2D:  Kinetostatic  Synthesis  of  Planar  Compliant  Mechanisms 

Objective:  In  this  task,  we  extend  the  capabilities  of  the  DAS  2D  mechanism  simulation  software 
that  allows  users  to  perform  kinematic  and  kinetostatic  synthesis.  Geometric  properties  and  com¬ 
pliances  of  the  flexible  members  are  synthesized  for  predefined  kinematic  and  static  requirements 
during  the  kinetostatic  synthesis.  Earlier,  we  had  developed  the  object  oriented  planar  mechanism 
analysis  tool  and  this  study  expands  the  capabilities  of  DAS  2D  software  to  include  interactive 
kinematic  synthesis.  We  also  present  a  solver  that  enables  users  to  access  and  use  main  analysis 
classes  of  the  software.  These  classes  can  be  used  for  custom  analysis  and  synthesis  cases  and  can 
be  integrated  to  other  software  programs. 

Methodology:  Kinematic  synthesis  is  a  systematic  attempt  to  calculate  the  design  variables  of 
a  mechanism  in  order  to  achieve  the  desired  function.  Static  analysis  should  be  combined  with 
the  kinematic  analysis  to  reveal  the  relationship  between  the  design  loads  and  the  accompanying 
deflection  of  a  compliant  mechanism.  This  relation  can  be  employed  in  kinetostatic  analysis  of 
compliant  mechanisms  for  obtaining  the  desired  relation  between  the  input  and  output  loads  and 
the  resulting  deflections. 
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There  are  a  number  of  software  programs  available  for  kinematic  synthesis  of  mechanisms. 
These  tools’  capabilities  are  mostly  limited  to  few  mechanism  types.  Linkages  [66]  is  capable 
of  planar  four  bar,  five  bar  and  six  bar  mechanism  synthesis  and  planar  mechanism  synthesis 
with  Sphinx  [67]  is  restricted  to  mechanisms  with  four  pin  joints.  Similarly,  kinematic  synthesis 
with  six  pin  joints  is  possible  [68].  Serial  chain  spatial  mechanism  synthesis  is  available  with 
Synthetica  [69]  and  likewise,  Spades  [70]  can  be  used  for  spatial  mechanisms  consisting  of  four 
cylindrical  joints.  Parallel  manipulators  with  RRR  and  RPR  legs  can  also  be  synthesized  [71]. 
SAM  [45]  is  a  commercial  mechanism  design  tool  and  is  capable  of  kinematic  synthesis  of  planar 
mechanisms.  Ch  Mechanism  Toolkit  [72]  enables  users  to  create  and  analyze  specific  mechanisms 
with  the  Ch  language.  Kinematic  synthesis  can  be  achieved  by  developing  case  by  case  synthesis 
codes.  Yu  et.  al.  [50]  also  developed  a  kinematic  synthesis  tool  that  can  synthesize  a  mechanism 
that  will  trace  the  input  coupler  curve.  However,  the  solver  is  limited  to  the  mechanisms  that  are 
stored  in  their  database. 

No  software  tool  is  available  for  synthesis  of  compliant  mechanisms,  although  the  fundamental 
theory  has  been  developed  for  some  time.  Topology  optimization  [73,74]  and  kinetostatic  synthesis 
[75]  are  the  two  principal  kinetostatic  synthesis  techniques  for  compliant  mechanisms.  Topology 
optimization,  which  is  more  suited  to  truss  type  mechanisms,  employs  finite  element  methods  to 
minimize  the  total  strain  energy  of  the  structure.  Kinetostatic  synthesis  methods  convert  the 
compliant  mechanism  to  rigid  body  mechanisms  via  pseudo-rigid-body  models  (PRBM)  [30]  and 
solve  the  synthesis  problem  using  the  resulting  rigid  body  mechanism. 

We  previously  developed  a  graphical  mechanism  simulation  tool  [65,76]  capable  of  unified  design 
and  analysis  of  planar  rigid  and  compliant  mechanisms.  As  a  next  step,  interactive  kinematic  and 
kinetostatic  synthesis  capabilities  were  developed.  The  software  is  developed  with  MATLAB  em¬ 
ploying  complete  object  oriented  approach.  Fig.  18  demonstrate  the  design,  analysis  and  synthesis 
classes  of  the  software.  The  graphical  user  interface  has  access  to  all  three  class  groups  and  provides 
interactive  design,  analysis  and  synthesis  of  mechanisms.  We  also  developed  a  solver  program  that 
allows  users  to  integrate  design  and  analysis  classes  to  other  MATLAB  applications.  These  classes 
allow  users  to  easily  create  and  analyze  any  planar  mechanism  with  a  few  line  of  codes.  Therefore, 
the  solver  broadens  the  capability  of  software  to  the  complex  cases  beyond  the  user  interface.  This 
paper  presents  the  interactive  synthesis  modules  and  the  solver. 

The  kinematic  synthesis  problem  is  defined  as  matching  the  coupler  curve  of  a  designed  mech¬ 
anism  with  a  desired  coupler  curve.  Type-P  [77]  and  Type-Z  [78]  are  two  similar  methods  for 
calculating  the  Fourier  descriptors  of  closed  or  open  curves.  Type-P  descriptors  are  calculated 
using  exact  positions  of  the  points  on  a  curve  and  on  the  other  hand,  Type-Z  descriptors  utilize 
the  slope  between  every  two  points  on  the  curve.  Type-P  descriptors  are  sensitive  to  scaling,  rota¬ 
tion  and  translation  resulting  in  difficulty  in  converging  to  the  optimal  solution.  A  point  and  the 
corresponding  slope  on  a  curve  are  defined  as: 


Pi  =  Xi  +  iyi  i=l,...,n  Si 


Pi+i  ~  Pi 

\\pi+i-pi\\ 


i=l,...,n-l 


where  n  is  the  number  of  points  on  the  curve. 

Type-P  descriptors  can  be  calculating  based  on  previously  defined  slope  information: 


F.D.k 


1 

n 


Tl—l 

Y  s-j  exp 

3= 0 


-  N  <  k  <  N 


where  N  determines  the  accuracy  of  the  curve  representation. 


(8) 


(9) 
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SOLVER 


Figure  18:  The  structure  of  the  DAS  2D  program. 


With  the  Type-P  descriptors,  the  kinematic  synthesis  problem  can  be  defined  as  an  constrained 
optimization  problem: 


/27V+1 


mm/=  (  \\RD-iM-RD-ihc)\ 


V  2=1 


subject  to  Optimization  Boundary 


(10) 


where  ^  are  the  link  lengths. 

Since  Eqs.(8-10)  are  not  sensitive  to  curve  translation,  rotation  and  scaling,  the  resulting  mech¬ 
anism  must  be  post-processed.  Rotation  and  scaling  can  lead  to  mechanism  moving  out  the  opti¬ 
mization  boundary  but  the  size  of  the  mechanism  will  not  be  altered.  However,  the  scaling  operation 
will  definitely  modify  the  size  of  the  optimized  mechanism.  The  scaling  operation  will  be  highly 
undesirable  if  a  very  small  part  of  the  coupler  curve  traces  the  target  open  curve  and  the  resulting 
mechanism  will  be  much  larger  than  the  optimization  boundary.  Due  to  this  undesirable  scaling  in 
open  curves,  the  slope  in  Eq.(8)  is  modified  as: 

Si  =  Pi+1  -  Pi  i=l,  n-1  (11) 

With  this  modification,  the  Type-P  transform  will  be  sensitive  to  scaling  operation  and  there¬ 
fore,  the  optimized  mechanism  does  not  need  to  be  scaled  after  optimization. 

A  kinetostatic  synthesis  case  arises  when  loads  acting  on  the  mechanism  are  known  and  a 
specific  deflection  is  desired.  In  this  case,  the  role  of  the  tool  is  determining  the  cross-section  (I) 
and  material  (E)  properties  of  the  flexible  members  in  the  mechanism.  In  this  case,  overall  design 
of  the  mechanism  is  known  but  the  flexural  stiffnesses  (El)  of  the  compliant  members  are  unknown. 
Therefore,  the  number  of  unknowns  are  equal  to  the  number  of  compliant  members.  The  synthesis 
problem  can  be  converted  to  the  following  optimization  problem: 

min  /  =  (Ades  -  Acur{EI ))  (12) 
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Bistable  mechanisms  are  special  kind  of  mechanisms  that  have  two  equilibrium  points.  The 
mechanism  will  stay  at  these  points  without  requiring  any  external  force.  Compliant  bistable 
mechanisms  acquire  two  stable  states  by  storing  energy  at  their  compliant  segments.  Critical 
load  is  defined  as  the  maximum  external  load  (force  or  moment)  that  is  required  to  actuate  the 
mechanism  from  one  stable  to  the  other  stable  position.  We  can  define  two  different  synthesis  cases 
for  bistable  mechanisms.  First  synthesis  mode  (Optimization  Variables:EA,  Eq.(13)  is  for  finding 
the  required  stiffness(es)  for  the  compliant  beam(s)  that  will  result  in  specified  critical  load. 

min  /  =  (C.L.des  ~  C.L.cur(EI ))  (13) 

El 

Alternatively,  in  the  second  synthesis  (Optimization  Variables:  node(a^,^),  Eq.(14)  desired 
instable  and  stable  positions  of  the  mechanism  can  be  determined  and  the  link  lengths  can  be 
found  via  optimization.  First  synthesis  mode  is  a  kinetostatic  synthesis  problem,  whereas  the 
second  synthesis  mode  is  a  kinematic  synthesis  problem. 


min  /  =  ( S.P.des  -  S.P.cur(node(x,y )))  (14) 

node(x,y) 


Constant  force  mechanisms  are  the  mechanisms  that  maintain  a  constant  output  force  during 
a  deflection  range.  Compliant  members  or  compliant  joints  can  also  be  used  in  a  constant  force 
mechanism.  For  a  compliant  mechanism,  the  link  length,  torsional  spring  stiffness  and  flexural 
stiffness  ratios  must  be  established  in  a  way  that  the  load  output  of  the  mechanism  will  be  almost 
constant.  The  optimization  (Eq.(15)  problem  can  be  stated  as: 


min  f  — 
Ri 


f  Fmax(Ri)\ 
\Fmin(Ri )  ) 


subject  to  Ibi  <  Ri  <  ub{ 


(15) 


Summary  of  Main  Results  The  capabilities  of  the  DAS  2D  early  design  tool  has  been  extended 
to  compliant  mechanism  synthesis.  There  are  four  (Fig.  19)  synthesis  modules  that  are  capable 
of  kinematic  and  kinetostatic  synthesis  of  compliant  mechanisms.  Kinematic  synthesis  module 
will  match  the  current  coupler  curve  of  a  mechanism  to  a  desired  curve.  Kinetostatic  synthesis 
module  will  find  the  required  flexural  stiffnesses  of  compliant  members  that  will  result  in  the 
prescribed  deflection.  Bistable  mechanism  synthesis  module  modifies  the  bistable  behavior  of  a 
bistable  mechanism.  Constant  force  mechanism  will  alter  the  design  of  a  mechanism  in  order  to 
have  a  constant  force  output. 

In  addition  to  the  synthesis  modules,  a  solver  that  can  be  integrated  to  MATLAB  is  developed. 
This  solver  enables  developers  to  access  the  design  and  analysis  classes  directly  in  MATLAB  and 
by  this  way,  the  developers  can  utilize  the  classes  in  their  codes  to  develop  custom  analysis  or 
synthesis  methods  for  compliant  mechanisms. 
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c)Bistable  Mechanism  Synthesis 


d)Constant  Mechanism  Synthesis 


Figure  19:  Four  different  synthesis  modules  are  shown. 
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Table  2:  Type  Synthesis  of  Flapping  Mechanisms  I 
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Table  3:  Type  Synthesis  of  Flapping  Mechanisms  II 


Research 

Group 


Project 


Kinematic  Diagram 


Adjacency  Matrix 


Characteristic 
Polynomial  Coeff, 


DelFly  I 


-110  1 
10  10 

0  10  1 
10  10, 


Delft 

University 


DelFly  II  & 
DelFly  Micro 


Y 


A 


-1  1  0  1 
10  10 
0  10  1 

1  0  10 


New 
DelFly  (I 


-I  1  0  1 
10  10 
0  10  1 
10  10, 


Cj-1 

Ci=l 

c2=-4 

Ci=-2 


c4=l 

c2=-4 

Ci=-2 


c4=l 
C3— 1 
C2=-4 
ct--2 


HMF 


'tLlT' 


-110  2 
10  2  0 
0  2  0  2 

2  0  2  0 


Harvard 

University 


PARITY 


w—f 


£ 


jT  • 


-1 

1 

0 

0 

2i 

1 

0 

2 

0 

0 

0 

2 

0 

2 

0 

0 

0 

2 

0 

2 

2 

0 

0 

2 

oJ 

t4-l 

c.=l 

c2--13 

Ci=*8 

c0=4 


c5=T 

C4=T 

c3=17 

C*=12 

c^-SG 

Cq-16 


University 

of 

Delaware 


FWMAV 


flapping  6 /  M  rotfctr 

spring  U ... 

■  connecting  rod 


-1  1  0  2 

10  10 
0  10  1 

2  0  10 


c4=l 

C3=l 

C2=-7 

Ci=-2 

c0=l 


Wing  Supports 

/ 


Small  Bird 


Compliant 

Frame 

Rocker 


Crank 


-110  0  3 

10  10  0 
0  10  10 
0  0  10  1 
3  0  0  1  1 


University 

of 

Maryland 


PmaifliK-jeca^ 


Jumho  Bird 


Ki|iidl  TicnIv 
Itvdme 

IdllH' 


-1 

1 

0 

1 

0 

2 


10  10 
0  10  0 
10  10 
0  10  2 
0  0  2  0 
0  0  0  2 


C'mnik 


C5=l 

C4=0 

c3=14 

c,=0 

c^-23 

cd=6 


c,=l 

c5=l 

c4=-1G 

ciF-io 

c2=32 

Ci~12 

c0=T6 
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Table  4:  Optimization  results 


Objective 

Constraint 

Optimized  variables 

Optimal  value 

%  Change 

h 

u 

w 

a 

A0g 

Max  L 

P  <  0.745 W 

6mm 

10mm 

1mm 

0.2527rad 

lrad 

48.95g 

39.36% 

Min  P 

L  >  35.12$ 

3.6mm 

10mm 

1mm 

0.2527rad 

lrad 

0.4044W 

-45.7% 

Min  P 

L  >  30g 

2.2mm 

10mm 

1mm 

0.2527rad 

lrad 

0.2151W 

-71.12% 

Max  L 

P  <  0.8W 

6mm 

9.7mm 

1mm 

0.2527rad 

lrad 

50.93g 

44.99% 

Max  L 

P  <  0.9W 

6mm 

9.2mm 

1mm 

0.2527rad 

lrad 

54.72g 

56.38% 

Table  5:  PRB  models  with  2  segments 


PRB 

model 

No. 

Schematic  DOF  of 

Symbolic  D 

Optimal  D 

f(fi) 

param 

1R  [30] 

1 

2 

oo 

[h 

OO 

OO 

7 

1  ~l\ 

PR 

2 

3 

OO 

ke 

kex 

OO 

7 

1  -7 

RP 

2 

3 

OO 

_ke 

OO 

k 

^ex 

l - 1 

O' 

^  1 
t— 1 

PRP  [39] 

3 

4 

OO 

ke 

kexl 

kex2 

7 

1  -7. 

oo 

oo 

0.451' 

0.934 

oo 

0.549 

oo 

0.012 

0.441 

0.942 

oo 

0.551 

oo 

oo 

0.465 

0.931 

0.007 

0.535 

oo 

1.361 

0.450 

0.934 

2.341 

0.550 

0.0396 

0.0366 

0.0321 

0.0395 
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Table  6:  PRB  models  with  3  segments 


PRB 

model 


Schematic  DOF 


No. 

of 


Symbolic  D 


Optimal  ft 


m 


/ 

OO 

oo 

7i" 

oo 

OO 

0.138' 

2R  [3?] 

2 

4 

ke  2 

oo 

72 

2.371 

oo 

0.640 

0.0036 

ke  3 

oo 

73. 

1.700 

oo 

0.222 

oo 

oo 

7i" 

oo 

oo 

0.139' 

RPR  [40] 


PRPR 


RPRP 


PRRP 


PRPRP 


k'0‘2  kex 

ke  3  oo 


72 
73. 

oo  kex  i  71 

ke  2  &ea;2  72 

OO  73 

oo  71 
k$2  kex  2  72 

&03  kex3  73 
OO  ^ea;l  7l 
ke2  OO  72 
&03  &ea;3  73. 

oo  kex\  7i 

&02  &ea;2  72 

^03  kex3  73 


.^03 

OO 


2.369 

1.701 

oo 

2.369 

1.700 

oo 

2.256 

1.759 

oo 

2.368 

1.701 

oo 

2.350 

1.711 


3.052 

oo 

2.987 

1.778 

oo 

oo 

0.099 

1.082 

3.096 

oo 

3.265 

1.716 

3.640 

0.679 


0.639 

0.222 

0.139: 

0.639 

0.222 

0.160: 

0.617 

0.222 

0.139: 

0.639 

0.222 

0.142: 

0.637 

0.221 


0.0036 


0.0035 


0.0031 


0.0035 


0.0034 


Table  7 :  PRB  models  with  4  segments 


PRB 

model 


Schematic 


3R  [32] 


RPRR 


RRPR 


RPRPR 


DOF 


3 


4 


4 


5 


No. 

of 

par  am 


6 


7 


7 


8 


Symbolic  O 


OO 

OO 

7i" 

ke  2 

OO 

72 

ke  3 

oo 

73 

ke  4 

oo 

74. 

oo 

oo 

7i 

ke  2 

kex2 

72 

ke  3 

OO 

73 

ke  4 

oo 

74 

oo 

oo 

7i 

ke  2 

oo 

72 

ke  3 

kex3 

73 

ke  4 

OO 

74 

oo 

oo 

7i 

ke  2 

kex2 

72 

ke  3 

kex3 

73 

ke  4 

OO 

74 

Optimal  D  f(D) 


oo 

3.491 

2.856 

2.716 

oo 

3.479 

2.881 

2.700 

oo 

3.130 

2.962 

2.880 

oo 

3.359 

2.828 

2.827 


0.095' 
0.392 
0.351 
0.161 
0.096' 
0.390 
0.352 
0.162 
0.1 14: 
0.409 
0.324 
0.153 
0.108: 
0.385 
0.353 
0.154 


0.0014 


0.0013 


0.0013 


0.0012 


oo 

oo 

oo 

oo 

oo 

3.395 

oo 

oo 

oo 

oo 

0.725 

oo 

oo 

3.188 

0.317 

oo 
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